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Theoretical studies within the chiral unitary approach, and recent experiments, have provided 
evidence of the existence of two isoscalar states in the region of the A(1405). In this paper we use 
the same chiral approach to generate energy levels in a finite box. In a second step, assuming that 
these energies correspond to lattice QCD results, we devise the best strategy of analysis to obtain 
the two states in the infinite volume case, with sufficient precision to distinguish them. We find 
out that using energy levels obtained with asymmetric boxes and/or with a moving frame, with 
reasonable errors in the energies, one has a successful scheme to get the two A(1405) poles. 



I. INTRODUCTION 

The history of the A(1405) as a composite state of 
meson baryon, dynamically generated from the meson 
baryon interaction, is rather long, starting from the works 
of Refs. [E Hi- Early works using the cloudy bag model 
also reached similar conclusions Q ■ The advent of chiral 
unitary theory, combining chiral dynamics and unitarity 
in coupled channels, brought new light into this issue and 
the A(1405) was one of the cleanest examples of states 
dynamically generated within this approach [^j-ffil . Hints 
that there could be two states rather than one had also 
been reported using the cloudy bag model Q and the chi- 
ral unitary approach ^ qualitative step forward 
was done in Ref. pT| , where two different versions of the 
approach were used, the two poles remained, and their 
origin was investigated. It was found that in an SU(3) 
symmetric theory there were two degenerate octets and 
a singlet of dynamically generated resonances, but with 
the breaking of SU(3) the degeneracy was removed, one 
octet with isospin 1 = moved to become the A(1670) 
and the other one moved close to the singlet, producing 
two poles close by in the region of the A(1405). One of 
the poles appears at energies around 1420 MeV, couples 
mostly to KN and has a small width of around 30 MeV. 
The other pole is around 1395 MeV, couples mostly to 
ttS and is much wider, around 120 or 250 MeV depend- 
ing on the model. After the work of Ref. [ll|, all further 
works on the chiral unitary approach have corroborated 
the two poles, with remarkable agreement for the pole at 
higher e nerg y and larger variations for the pole at lower 
energies ^^U^. 

Suggestions of experiments to confirm this finding were 
made, and it was shown that one should not expect to 
see two peaks in the cross sections, but rather differ- 
ent shapes in different reactions. In this sense, a sug- 
gestion was made to look for the A(1405) peak in the 
K^p — > 77rS reaction [lO], where the 7 would be radi- 
ated from the initial state, making the K~p system to 



lose energy and go below threshold and then excite the 
high energy state of the A(1405), to which it couples 
most strongly. This reaction was not made, although 
it is planned for JPARC [2l[, but a similar one, where 
the photon was substituted by a pion, was implemented 
in Ref. [l^l studying the K^p — > tt^tt^TP reaction at 
Pk = 514 MeV/c - 750 MeV/c. A neat and narrow 
peak was seen at ^/s = 1420 MeV, which was analyzed in 
Ref. [2^ and interpreted in terms of the high energy pole 
of the A(1405). More recently it was noticed that old data 
on the K~d — >■ ttS n reaction from Ref. [2^] produced 
a peak in the ttS spectrum around y/s = 1420 MeV, 
with also a small width. These data were well repro- 
duced in Ref. (25j within the chiral unitary approach and 
multiple scattering, and once again it was shown that 
it gave support to the existence of the second pole of 
the A(1405). It was shown there that the reaction pro- 
ceeded with kaons in flight but not for stopped kaons, 
because the background from single scattering was too 
large in this latter case, obscuring the signal of the res- 
onance that stems from double scattering. Even then, 
it was shown in Ref. [26J that kaons from the DAFNE 
facility, coming from the decay of the (j>, would also be 
suited to search for this resonance if neutrons were mea- 
sured in coincidence in order to reduce the background. 
The search for reactions where the A(1405) is produced 
has continued, showing that, as predicted, different reac- 
tions have different shapes. In this sense there have been 
recent photoproduction experiments [27l . [28j and proton 
induced experiments [l^llo] where the shapes are indeed 
different and the peaks appear at lower energies, around 
1405 MeV, as the nominal mass. There are also theoret- 
ical studies for these reactions where the peaks appear 
around these energies, and the larger contribution of the 
lower energy state that couples mostly to ttE is mostly 
responsible for it [31L35]. 

In as much as chiral dynamics is a good representa- 
tion of QCD at low energies, the predictions of the chiral 
unitary approach on the A(1405) stand on firm ground. 
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Yet, it would also be very interesting to have these pre- 
dictions confirmed with lattice QCD simulations. In this 
sense, the determination of hadron spectra is one of the 
challenging tasks of Lattice QCD and many efforts are 
being devoted to this problem [36l - l53 |. some of them in 
particular to the search of the A(1405) [55l - l60l |. A review 
on the Lambda(1405) and attempts to see it from differ- 
ent points of view is given in Ref. [6l| . In some works 
the "avoided level crossing" is usually taken as a signal 
of a resonance, but this criteria has been shown insuffi- 
cient for resonances with a large width [g^ - Ig^ . Some- 
times, the lattice spectra at finite volumes is directly as- 
sociated to the energies of the states in infinite volume 
invoking a weak volume dependence of the results, as 
done recently searching for the A(1405) resonance [6a]- 
A more accurate method consists on the use of Liischer's 
approach, for resonances with one decay channel. The 
method allows to reproduce the phase shifts for the de- 
cay channel from the discrete energy levels in the box 
(60 . [13] ■ This method has been recently simplified and 
improved in Ref. (63 ] by keeping the full relativistic two 
body propagator (Liischer's approach keeps the imagi- 
nary part of this propagator exactly but makes approx- 
imations on the real part) and extending the method to 
two or more coupled channels. The method has also been 
applied in Ref. [68[ to obtain finite volume results from 
the Jiilich model for meson baryon interaction, includ- 
ing spectra for the A(1405) with finite volume, and in 
Ref- 169], to study the interaction of the DK and rjDs 
system where the Ds*o(2317) resonance is dynami cally 
generated from the interaction of these particles [70l - [73| . 
The case of the k resonance in the Ktt channel is also 
addressed along the lines of Ref. [gl] in Ref. [tI] . 

In the work of Ref. [gj , the inverse problem of getting 
phase shifts and resonances from lattice results using two 
channels was addressed, paying special attention to the 
evaluation of errors and the precision needed on the lat- 
tice results to obtain phase shifts and resonance proper- 
ties with a desired accuracy. Further work along these 
lines is done in Ref. [t^ . The main problem encountered 
is that the levels obtained from the box of a certain size 
range do not cover all the desired energy region that one 
would like to investigate. Several suggestions are given in 
order to produce extra levels, like using twisted bound- 
ary conditions or asymmetric boxes [6^ . These are, how- 
ever, not free of problems since it is unclear whether a full 
twisting can be done in actual QCD simulations including 
sea quarks, and the asymmetric boxes have the problem 
of the possible mixing of different partial waves. Another 
alternative is to evaluate levels for a system in a moving 
frame as done in Ref. [53], but this also poses problems 
of mixing in principle. The generalization of Liischer's 
approach to the moving frame is done in Rcfs. 75 79], 
and it provides a convenient framework for lattice calcu- 
lations since new levels can be obtained without enlarging 
the size of the box, with an economy in computational 
time. It is then quite convenient to carry out simulations 
using effective theories in a finite volume, preparing the 



grounds for future lattice calculations, trying to find an 
optimal strategy on which configurations to evaluate in 
order to obtain the desired observables in the infinite vol- 
ume case. 

The case of extracting the A(1405) parameters is spe- 
cially challenging, particularly because two resonance 
must be found which are not too far from each other, 
which means that extra precision will be demanded to 
the lattice results. Furthermore, the two poles are not to 
be seen in the ttS phase shifts, since, as mentioned before, 
different amplitudes give different weight to the two poles 
and the ttE phase shifts provide insufficient information. 
The other reason is that the chiral unitary approach tells 
us that the two states couple strongly to KN and ttS, so 
the use of the two channels in the analysis is mandatory 
and the use of one channel like in Luscher approach is 
bound to produce incorrect results. In view of this we 
face the problem using the two channels explicitly in the 
analysis and produce amplitudes in the coupled chan- 
nels from where we can extract the pole positions in the 
complex plane by means of an analytical continuation 
of these amplitudes. Even then, the problem is subtle 
because using standard periodical boundary conditions, 
and a wide range of lattice volumes, there is a gap of 
energies in the levels of the box precisely for the energies 
where one finds the poles. Because of this problem one is 
then forced to use either asymmetric boxes or discretiza- 
tion in the moving frame in order to find eigenvalues of 
the box in the desired region. In the present paper we 
face all these problems and come out with some strate- 
gies that we find better suited to determine the position 
of the two A(1405) poles. 



II. FORMALISM 

In the chiral unitary approach the scattering matrix in 
coupled channels is given by the Bethe-Salpeter equation 
in its factorized form 



T = [1 - VG]-W = [V-'^ - G]-\ 



(1) 



where V is the matrix for the transition potentials be- 
tween the channels and G is a diagonal matrix with the 
i^^ element, Gi, given by the loop function of two prop- 
agators, a pseudoscalar meson and a baryon, which is 
defined as 



G, = i2M, 



d^p 1 
(27r)4 (P - p)2 _ m2 + j£ ^ 



1 



(2) 

where rrii and Mi are the masses of the meson and the 
baryon, respectively, and P the four-momentum of the 
global meson-baryon system. 

The loop function in Eq. ^ needs to be regularized 
and this can be accomplished either with dimensional 
regularization or with a three-momentum cutoff. The 
equivalence of both methods was shown in Refs. [8l[8l|. 
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In dimensional regularization the integral of Eg. ^ is 
evaluated and gives for meson-baryon systems [Sj, [82] 



Gi{s,mi, Mi) = 



(47r)2 



TTl- 

a,(^) + log^ 



+ 



2s 



-log^ 



log(s-(M2-TO2) + 2ViQ.(Vi)) 

+ log [s + (Aff - m^) + 2ViQz(V^)) 
- log (-S + (Mf - mj) + 2^sQ,{V~s)) 

_ log (-S - (Mf - m^) + 2^sQ,{^s)) 



(3) 



where s — E^, with i? the energy of the system in the 
center of mass frame, Qi the on shell momentum of the 
particles in the channel i, fi a regularization scale and 
ai(fi) a subtraction constant (note that there is only one 
degree of freedom, not two independent parameters). 

In other works one uses regularization with a cutoff in 
three momentum, once the p'^ integration is analytically 
performed [83], and one gets 



IpI <Pmax 

Wl,2(p) = 



2M,: 



UJi{p) + UJ2{P) 



(2^)3 2LUiip)LU2ip) - (Luiip) + LJ2{p))^ + ie' 



(4) 



with mi, m2 corresponding to mi and Mi of Eq. ([2|). 

When one wants to obtain the energy levels in the fi- 
nite box, instead of integrating over the energy states of 
the continuum, with p being a continuous variable as in 
Eq. ((?]), one must sum over the discrete momenta allowed 
in a finite box of side L with periodic boundary condi- 
tions. We then have to replace G by G = diag (Gi, G2) 
(in two channels), where 



1 ^i(p) -f- UJ2{'P) 

L3 ^ 2uji{p) uj2{p) E^ - {uji{p) + L02{p)Y ■ 



^ 27r ^ _^ . 
p = — n, n G Z 
L 



(5) 



This is the procedure followed in Ref. [6J]. The 
eigenenergies of the box correspond to energies that pro- 
duce poles in the T matrix, Eq. ([T]) , which correspond to 
zeros of the determinant of 1 — VG, 

det(l - V^G) = 0. (6) 

For the case of two coupled channels Eq. ^ can be writ- 



ten as 



det(l - VG) 



^ 1 - ViiGi 

+ (VilV^22 - 

: 0. 



V22G2 

r2 - 



Vi2)GiG2 



(7) 



The problem of the KN interaction with its coupled 
channels and the A(1405) was addressed in Ref. ^Sj us- 
ing the cut off method, but more recently it has been 
addressed using dimensional regularization [1, [s^] • For 
this reason we will also use the dimensional regulariza- 
tion method for the finite box, which was developed in 
Ref. [g^ . The change to be made is also very simple, the 
G function of dimensional regularization of Eq. ^ has 
to be substituted by 



G{E) = G"{E)+ lim 



^3 E 



L3 



_fp_ 

(27r)3 



Hp) 



D<p„ 



where I{p) is given by 



Iip) = 



2M, 



LOlip) + U!2{P) 



2uJi{p) uj2ip) E^ - {uji{p) +uj2ip))^ +«e' 



(8) 



(9) 



We will also consider the case where the meson-baryon 
system moves with a fourmomentum P = (P°, P) in the 
box. In this case we still have to define the integrals and 
the sums in the CM frame, where Pmax is defined, but 
the momenta of the two particles must be discretized in 
the box, where the system moves with momentum P. We 
follow the approach of Refs. [t^, ^8d\ and use the boost 
transformation from the moving frame, with the variables 
Pi, 2, to the CM frame with the variables Pi2 



Pl,2 =Pl,2 



Mi \ pi,2 ■ P P*^2 

pO J \p\2 PO 



P. 



(10) 



where Mf = P"^ — P"^ ^ P^, the subindexes 1,2, rep- 
resent the meson, baryon particles and ^^"2 the CM 
energies of the particles given by 



Pl,2 



Mf + ml 2 - m\ 



2M 



J/ 



Then we must do the substitution in Eq. 
evaluation of the energies in the box. 



(11) 
for the 
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|p*|<Pmax 

1 Ml ^ 



L3 



27r 



p = —r'^i € 



(12) 
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FIG. 1. Energy levels in a symmetric box of side length L. 

with p* given in terms of pi by means of Eq. pup . 

Since pi andp2 — P^Pi must both satisfy the periodic 
boundary conditions, this forces P to be also discretized 
and thus we can only use values of P such that 

- 27r ^ ^ . 
P = —N, N ^1? . (13) 



III. RESULTS 

A. Energy levels in the box 

In this section we show the energy levels obtained from 
the solution of Eq. ([5]) as a function of the side length of 
the box, L, and for different physical cases: using peri- 
odical boundary conditions in a: (1) symmetric box, (2) 
asymmetric box and (3) symmetric box but in a moving 
frame, i.e., with non-zero value for the total center of 
mass momentum P (Eq. (|13|1 ). 

1. Periodical boundary conditions in a symmetric box 

In Fig. [U we show the first six energy levels related to 
the system formed by the coupled channels KN ^ ttE, r]K 
and K'E., which generate a double pole structure for the 
A(1405) and a pole for the A(1670) These levels 

are obtained by solving Eq. ^ using the chiral model of 
Ref. [m and imposing periodical boundary conditions in 
a symmetric box of side length L (measured in units of 

As it can be seen in Fig. [TJ the gap between the levels 
0, 1 and specially between levels 1 and 2 is considerable, 
giving rise to the presence of only two levels in the en- 
ergy region of interest, i.e., the energy range in which the 
two poles of the A(1405) are found (1390 - 1430 MeV). 
This fact shows the difficulty that one can face to extract 



FIG. 2. Energy levels in an asymmetric box of side length 
Lx = Ly = L and Lz — zL, with z — 0.5 — 2.5 in steps of 0.5. 



information about the poles of the A(1405) in an infinite 
volume considering these energy levels as reference. 



2. Periodical boundary conditions in an asymmetric box 

To see if we can obtain more energy levels in the region 
of the A(1405), it is also possible to solve Eq. (O but in 
an asymmetric box. To do this we just need to substitute 
by Lz and the momentum p of Eq. by p = 
{2Tr)(rix/Lx,ny/Ly,nz/Lz)- In Fig. [2 we show the first 
three energy levels determined in a box of side lengths 
Lx = Ly = L and Lz = zL, and we vary z between 0.5L 
and 2.5L. In this way, we get more energy levels in the 
region of interest, which can provide different information 
about the system and the poles of the A(1405). 



3. Periodical boundary conditions in a moving frame 

Another method to try to get more energy levels 
around the pole positions of the A(1405) and thus, differ- 
ent information about the dynamics of the system under 
consideration, consist in imposing periodical boundary 
conditions in a symmetric box of side length L but con- 
sidering the system in a moving frame, i.e., with non zero 
center of mass momentum P. In Fig. [3] we show the re- 
sults found in this case for the first three levels obtained 
and for different values of the vector N (see Eq. (|13p ). As 
it can be seen, the use of different values of P gives rise to 
a splitting of the levels. In particular, the sphtting of level 
1 is precisely in the energy region of interest, 1390— 1450 
MeV. This is different to the case of the asymmetric box, 
where level 2 is required in order to have energy levels 
around 1420-1450 MeV, as it can be seen in Fig.O 
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FIG. 3. Energy levels in a symmetric box of side length L 
with the system having a center of mass momentum given by 
Eq. dSl). 



B. The inverse problem: getting the A(1405) poles 
from the energy levels of the box 

In the following we refer to the problem of determin- 
ing the pole positions of the A(1405) in the infinite vol- 
ume using the energy levels shown in Figs. (TJ [51 [3] as if 
they were provided to us by a lattice calculation. In our 
formalism, we can simulate lattice-like data considering 
points related to the energy levels of Figs. [TJ [21 [3] and 
assigning to them a typical error of ± 10 MeV. We call 
the data generated in this form "synthetic" lattice data 
and the problem of getting the poles of the A(1405) from 
these data points "the inverse problem". 

To solve the inverse problem we consider a poten- 
tial with the same energy dependence than the chiral 
potential used to generate the energy levels shown in 
Figs. [U [21 [21 In its non-relativistic version, this potential 
is given by ^GJ 



(14) 



with Cij coefficients depending on the channel consid- 
ered, / the pion decay constant and Ei {Ej) the center 
of mass energy of the meson in the initial (final) state. 
Using that for a particular channel / 



E, - 



E' 



Mf 



2E 



(15) 



with mi and Mi the masses of the meson and baryon 
which constitute the channel /, respectively, we can write 
Eq. (HID as 



-^{E+^hl+^^'iM! + M-)]\ (16) 



Choosing a region of energies around a certain value of 
E, Eq, the inverse function of E can be expanded as a 
function oi E — Eq to a good extend. Particularizing Eq 
to the value given by the sum of the kaon and nucleon 
masses, i.e., niK + M^, we can write the potential in 
Eq. 119 as 



bij[E - (rriK 



Mn)]. 



(17) 



The value of the coefficients a.^ and bij can be obtained 
comparing Eq. P?|) with Eq. and substituting 1/E 
by its Taylor expansion around Eq = rriK + Mn ■ 

To solve the inverse problem, we use the energy levels 
obtained from Eq. © with the potential of Eq. pTj) but 
treating a.y and hij as parameters which are determined 
by fitting the corresponding solutions for the energy lev- 
els to the "synthetic" lattice data considered. Since this 
potential has the same energy dependence than the chi- 
ral potential, the best fit we can perform will have as 
minimum value for the the result Xmin ~ 0- How- 
ever, other possible potentials, giving rise to solutions 
compatible with the error assumed in the data points, 
can be also found as an answer for the inverse problem. 
These solutions can be obtained by generating random 
numbers for the parameters and bij close to those of 
the minimum such that ^ Xmin + 1- 

It is important to notice that the loop function G, used 
in Eq. ([S]), needs to be regularized and, thus, depends on 
a cut-off or a subtraction constant. Consequently, so do 
the fitted parameters, but the T matrix obtained from 
Eq. IT]) and the observables related to it should be in- 
dependent of this regularization parameter. This means 
that the inverse method cannot depend on the cut-off 
or subtraction constant assumed in the evaluation of the 
G function. For the case of one channel, it is possible 
to show analytically this independence in the choice of 
the cut-off or subtraction constant [13, , but if more 
channels are involved it can only be seen numerically by 
changing the cut-off or subtraction constant in a reason- 
able physical range [H, [1^ . 

In the next sections we show the results found for the 
inverse problem. To accomplish this we have considered 
different sets of points extracted from the energy levels 
shown in Figs. [1] [21 [21 and fitted them from the solution 
that Eq. ^ produces with the potential of Eq. pT]) . To 
solve Eq. ^ we have taken into account two coupled 
channels, ttS (which we named channel 1) and KN (or 
channel 2), which are the most relevant channels to de- 
scribe the properties of the A(1405). This implies, as it 
can be seen in Eq. ([7]), that we have to determine three 
potentials, Vii, V12 (V21 — V12) and V22 or equivalently 
6 parameters an, ai2, a22, ^ii, bi2 and b22- Once the pa- 
rameters and, thus, the potentials, are known, we can use 
them to solve Eq. ([T]) and determine the pole positions 
of the A(1405) in an infinite volume. 
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1.5 2.0 2.5 3.0 

L(m-^ ) 

FIG. 4. First two energy levels as function of the box side 
length L, reconstructed from fits to the "synthetic" data of 
Fig. [T] in a range of L between 1.5 m^^ and 3.34 m^^ using 
the potential of Eq. p7p . The band corresponds to different 
choices of parameters within errors. 



1. Periodical boundary conditions in a symmetric box 

In Fig. 3] we show the resuhs of the energy levels recon- 
structed from the best fits to the "synthetic" lattice data 
considered from Fig. [TJ These data consist in 10 points 
for levels and 1 obtained in a symmetric box of side 
length L, varying L in the range 1.5 to 3.34 m^^, 
assigning an error of ± 10 MeV to the eigenenergies of 
the box (from now on, we will always assume an error 
of ± 10 for the different "synthetic" data that we will 
use). The shadowed band in the figure corresponds to 
the random choices of parameters satisfying the condi- 
tion ^ Xmin + 1- Using the potentials obtained from 
the fit and the loop function G in infinite volume, we 
can solve Eq. ([1]) and calculate the two-body T matrix 
in the unphysical sheet, which allows us to determine the 
pole position of the A(1405) associated to the band of 
solutions shown in Fig. 01 As a result we get a double 
pole structure for the A(1405), with one pole in the re- 
gion 1385-1433 MeV and half width between 93-137 MeV 
(which we call pole 1) and another one in the energy re- 
gion 1416-1427 MeV and half width in the range 11-20 
MeV (which we call pole 2). If we compare these results 
with the ones of the chiral model jjjj], 1390-i66 MeV and 
1426-il6 MeV, respectively, we find a big dispersion in 
the determination of the real part of the first pole of the 
A(1405). This shows that the information which one can 
extract from the "synthetic" data considered in Fig. 2] is 
not sufficient to determine with more precision the poles 
associated with the A(1405). 

A way to delimit the poles of the A(1405) with more 
precision from lattice data could consist in going to 
higher volumes, since for large volumes the results in the 
box should be very close to those of an infinite volume. 
With this idea in mind, we can generate "synthetic" data 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
L(m"^ ) 

FIG. 5. Same than in Fig. [4] but for a range of L between 
1.5 m^^ and 4.93 m~^. 

points for the levels and 1 of Fig.[Tl but in a larger range 
of L than the one considered in Fig. [31 The data points, 
as well as the results from the fits, are shown in Fig. [51 
Similarly, if we use now the potentials associated with the 
band of solutions shown in Fig.[51to solve Eq. ([T]) and cal- 
culate the T matrix in the unphysical sheet, we get again 
two poles in the complex energy plane associated with 
the A(1405): one in the region 1390-1433 MeV, with half 
width between 70-100 MeV, and other at 1410-1421 MeV 
with half width 17-30 MeV. Comparing them with the 
previous results, we find that the consideration of a big- 
ger box has improved slightly the width associated with 
the first pole of the A(1405), however, we continue having 
a similar energy dispersion for the real part of the pole. 

We could also try using different levels than the ones 
employed in Figs. [H and [51 to see if we can get more re- 
liable information from them. In Fig. [SI we consider a 
"synthetic" data obtained from levels 1 and 2 of Fig. [2 
We have taken into account 5 points for level 1 in a range 
of L between 1.5 m^^ and 3.9 to~^ and 4 points for level 
2 for values of L inside 2m~^ to 3.9 m~^. This is be- 
cause for level 2 the points for values of L below 2m~^ 
are influenced by the 77A and K'E channels and, thus, it 
is not possible to make a fit to them considering only the 
ttS and KN channels, as we do. We can use now the po- 
tentials associated with the different fits shown in Fig. [SI 
to calculate the pole positions of the A(1405) in infinite 
volume by means of Eq. ([T]). In this case, we continue 
getting a double pole structure for the A(1405), but this 
time one pole is at (1375- 1430) -i(70- 85) MeV and the 
other one is at (1412 - 1427) - i(21 - 34) MeV. The posi- 
tion of the second pole remains basically the same as in 
the two previous cases. However, the use of "synthetic" 
points generated from levels 1 and 2 instead than from 
levels and 1 has restricted more the imaginary part 
of the first pole, although we continue getting a similar 
energy dispersion for the real part of the pole position. 

Finally, we could consider all the energy levels present 
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3.5 4.0 



L(m^ ) 



L(m;; ) 



FIG. 6. Fits to the levels 1 and 2 of Fig. [T] constructed from 
the potential of Eq. (fT7|) . 



FIG. 8. Fits to the "synthetic" data extracted from the energy 
levels 0, 1 and 2 of Fig.[2]in an asymmetric box of side lengths 
= Ly = L and L^, — zL. The data points considered are 




L(m;; ) 



FIG. 7. Fits to the levels 0, 1 and 2 of Fig. [T] constructed 
from the potential of Eq. (|17p . 



in Fig.[T]below 1600 MeV to generate data points to check 
if the consideration of more levels can restrict more the 
energy region at which the pole positions of the A(1405) 
are found. Following this idea, in Fig. [7] we consider a 
set of 14 points extracted from levels 0, 1 and 2 of Fig.[T] 
Similar to the previous results, the consideration of data 
points associated to three energy levels puts a restriction 
on the imaginary part of the first pole of the A(1405), 
which in this case is in the range 54-68 MeV (closer to 
the chiral solution, 66 MeV). However the dispersion on 
the real part continues basically equal, 1400-1428 MeV. 
For the second pole we get (1408-1425)-i(29-40) MeV. 

These results show that the information which can be 
extracted from "synthetic" data constructed from the en- 
ergy levels obtained in a symmetric box of volume is 
not enough to determine with precision the pole posi- 
tions of the A(1405), a fact which is basically related to 
the absence of energy levels, and thus information about 



generated from level for z - 
z = 2.0 and level 2 for z = 2.0. 



2.5, level 1 for z = 0.5 and 



the dynamics of the system, in the region between 1400- 
1500 MeV, as it can be seen in Fig. [T] 



2. Periodical boundary conditions in an asymmetric box 

We consider now the case of an asymmetric box of side 
lengths Lx — Ly — L and = zL to solve the inverse 
problem. In this case, we generate a set of 20 data points 
extracted from levels 0, 1 and 2 shown in Fig. [21 In 
particular, we use 5 points for level calculated with 
z — 2.5, 10 points for level 1 (5 for the case z = 0.5 and 
5 more for z = 2.0) and 5 points for level 2 obtained 
with z — 2.0. In this way we ensure the presence of some 
energy levels in the region of the A(1405), as it can be 
seen in Fig. [51 

The solution of the Bethe-Salpeter equation in an in- 
finite volume, Eq. ([T]), using the potentials related to 
the band of solutions plotted in Fig. [51 shows the pres- 
ence of a double pole structure for the A(1405) with 
pole positions at (1383 - 1407) - i{57 - 69) MeV and 
(1425 - 1434) - i{25 - 35) MeV. Thus, using this new 
set of data points, there is an improvement in the deter- 
mination of the first pole of the A(1405), which is now 
quite close to the chiral result (1390-i66 MeV). However, 
the second pole appears at higher energies as compared 
to the case of a symmetric box and sometimes is far from 
the chiral solution (1426 — il6 MeV), being even close to 
the KN threshold. 



3. Periodical boundary conditions in a moving frame 

We can also study the information which can be ex- 
tracted for the poles of the A(1405) using the levels ob- 
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FIG. 9. Fits to the "synthetic" data extracted from the en- 
ergy levels and 1 of Fig. [31 which correspond to the case 
of a symmetric box, but with the particles being in a moving 
frame. 



tained when we consider the system in a symmetric box, 
but in a moving frame, to generate "synthetic" lattice 
data. In this case, we consider levels and 1 of Fig. [3] 
determined for 5 different values of the center of mass 
momentum (the ones shown in the legend of Fig. [S]) and 
two points in each of these curves. In particular, we take 
points at L = 1.757 to~^ and L = 2.014 m~^, obtaining 
then 20 data points. The results are shown in Fig. O 
From the solution of the best fits, we can use the poten- 
tials obtained to solve Eq. ([1]), getting then two poles for 
the A(1405): one at (1388 - 1418) - i{59 - 77) MeV and 
other at (1412 - 1427) - i{16 - 34) MeV. 

In Fig. [TO] we show the results for the pole positions of 
the A(1405) obtained from the different data set consid- 
ered in this work. As it can be seen in Fig. [TOl out of the 
different data sets considered to solve the inverse prob- 
lem, the cases of an asymmetric box and of a symmetric 
box but in a moving frame seem to be more suited to get 
the two poles of the A(1405) with more precision. 



IV. CONCLUSIONS 

We have made a study of the KN interaction with 
its coupled channels in a finite box and found the levels 
obtained as a function of the box size. We have done this 
for standard periodic conditions and symmetric boxes, 
for asymmetric boxes and for symmetric boxes but with 
the particles in a moving frame. The aim of the work 
has been to solve the inverse problem in which, assuming 
that the levels in the box would correspond to "QCD 
lattice results" we want to determine the pole positions 
in the complex plane for the two A(1405) states provided 
by the chiral unitary approach and supported by several 
experiments. 

We found that the problem is not trivial, and even 
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FIG. 10. Pole positions of the A(1405) reconstructed from 
the different set of "synthetic" data generated for the differ- 
ent cases considered in this work. The shaded symbols cor- 
responds to the positions obtained for the first pole of the 
A(1405), while the empty symbols are related to the second 
pole of the A(1405). 



the use of a large number of energies of the box corre- 
sponding to different levels and volumes with standard 
periodic conditions cannot provide the mass and width 
of the states with the accuracy of the chiral unitary ap- 
proach and present experiments. For this reason we in- 
vestigated other possible strategies and found that the 
use of asymmetric boxes and levels coming from the par- 
ticles in moving frames helped considerably to narrow 
down the uncertainties in the determination of the mass 
and width of these resonances. The choices of levels and 
energies made for this analysis should be a guiding tool 
for future QCD lattice evaluations, showing the num- 
ber of levels needed, the errors that should be demanded 
in the determination of the energies of the box and the 
type of asymmetric boxes or total momenta of the pair 
of particles in the moving frames. Having this informa- 
tion before hand is of tremendous value given the time 
consuming runs of actual QCD lattice runs. 
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